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INTRODUCTION 


Aircraft  structural  surfaces  often  serve  as  the  walls  of  integral  fuel  cells.  Projectiles 
penetrating  the  fuel  generate  intense  pressure  waves  capable  of  producing  catastrophic 
failure  of  the  ftiel  cell  walls.  This  phenomenon,  termed  hydraulic  ram,  can  lead  to  structural 
kill  of  the  aircraft.  Massive  fuel  >oes,  resulting  IVom  fuel  cell  failure,  can  also  lead  to  aircraft 
kill  by  on-board  fires  or  by  fuel  starvation^ 


'’Drag  forces  acting  to  decelerate  the  projectile  arc  the  source  of  pressure  waves  that 
propagate  into  the  fuel.  A  prominent  characteristic  of  the  fluid  flow  is  a  cavity  that  extends 
from  the  projectile  hack  to  the  impact  point.  The  cavity  grows  radially  away  from  the 
'  ^'-''  trajectory.  The  cavity  eventually  collapses, Imt  on  a  longer  time  scale  than  of  interest  here. 
The  fluid  pressure  field  is  highly  peaked  at  the  projectile.  The  fluid  velocity  field  is  high  at 
the  projectile  and  near  the  cavity. 


Pressure  waves  incident  upon  a  fuel  cell  wall  produce  reflected  waves  that  depend 
upon  the  motion  of  the  wall.  The  walls,  in  turn,  move  in  response  to  the  total  pressure 
loading.  The  coupling  between  the  fluid  pressure  field  and  wall  motion,  called  the 
fluid/structure  interaction,  is  of  critical  importance  with  dense,  condensed  fluids  like  fuel  or 
water.  In  contrast,  for  highly  compressible  fluids  like  gases,  the  fluid/structure  interaction 
can  be  neglected.  In  hydraulic  ram,  the  interaction  is  complicated  by  cavitation,  which  occurs 
in  response  to  negative  pressures  that  can  be  obtained  at  the  fluid/structure  interface. 

Motion  of  a  structure  can  be  readily  calculated  given  the  pressure  loading.  However, 
the  structural  response  to  hydraulic  ram  can  include  nonlinear  deflections  and  plastic 
deformations.  Fuel  cell  structure  can  be  either  metallic  or  fiber  composite.  The  construction 
may  include  elastomer  fuel  cell  liners.  Foam  or  air  gap  isolation  may  be  employed  between 
the  fuel  and  the  structure.  Of  interest  is  the  damage  that  occurs  during  hydraulic  ram:  walls 
may  fail  catastrophically,  joints  and  rivets  may  fail,  and  fiber  composites  may  delaminate  in 
response  to  bending. 

The  pressure  wave  generation  theory  had  its  roots  in  the  work  of  Bristow  and 
Lundeberg  (Reference  1)  and  Yurkovich  (Reference  2).  Yurkovich,  in  particular,  used  the 
idea  of  describing  the  fluid  motion  by  a  moving  point  source  in  a  potential  flow.  This  idea 
formed  the  basis  of  the  hydraulic  ram  pressure  theory  for  subsonic  tumbling  projectiles.  The 
theory  was  developed  by  Lundstrom  (Reference  3),  who  used  a  moving  line  of  sources  whose 
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strength  was  based  on  the  projectile  drag  characteristics.  The  theory  was  substantiated  by  a 
large  body  of  experimental  data  obtained  by  Lundstrom  and  Fung  (References  4  and  5). 

Predicting  structural  response  to  hydraulic  ram  pressure  toading  has  proven  to  be  a 
difficult  task.  Early  work  was  based  on  use  of  the  piston  theory  for  the  fluid/structure 
interaction.  This  theory  assumed  that  pressure  waves  reflected  from  fuel  tank  walls  were 
one-dimensional,  and  propagated  normally  away  from  the  wall  surface. 

Ball  (Reference  6)  modified  two  structural  analysis  computer  codes,  including  the 
piston  theory  and  hydraulic  ram  pressure  loading.  The  SATANS  code  was  limited  to  elastic 
motion  of  shells  of  revolution,  but  was  computationally  efficient. 

The  BR-1  code  (Reference  7)  was  a  more  general  finite  element  structural  code  that 
allowed  plastic  deformation  of  the  elements.  However,  this  code  was  very  computationally 
intensive,  taking  50  times  more  computer  time  than  the  SATANS  program.  The  BR-1  code 
version,  which  included  the  piston  theory,  was  called  the  BR-IHR  code.  Predictions  of  the  BR- 
1HR  and  modified  SATANS  codes  were  compared  with  experiments  in  Reference  8.  It  was 
found  that  experimental  strain  data  greatly  exceeded  the  predicted  values.  The  disagreement 
was  attributed  to  the  failure  of  the  piston  theory  for  hydraulic  ram  pressure  loading. 

The  Moir§  fringe  technique  was  used  by  Fry,  Newman,  and  Bless  (Reference  9)  to 
measure  the  deflection  of  entrance  panels  induced  by  hydraulic  ram.  The  data  were  analyzed 
by  Bless,  Fry,  and  Barber  (Reference  10)  using  the  BR-IHR  program,  with  results  similar  to 
those  of  Ball. 

Since  the  piston  theory  failed  to  describe  the  fluid/structural  interaction  in  hydraulic 
ram,  Lundstrom  (Reference  11)  proposed  the  variable  image  method.  The  fluid  potential  of 
pressure  waves  reflected  by  a  moving  wall  is  assumed  to  be  in  the  form  of  an  image  of  the 
incident  wave  potential  multiplied  by  an  unknown  function.  The  function  is  determined  by 
matching  the  fluid  velocity  with  the  wall  vel  icity  at  the  interface. 

At  the  same  time,  Ankeney  (Reference  12)  suggested  an  approximate  structural 
model  for  use  in  hydraulic  ram  analysis.  A  panel  deflection  profile  in  the  shape  of  a  cosine 
curve  was  assumed,  and  then  the  amplitude  was  calculated  using  the  Lagrange  equations  of 
motion. 


Lundstrom  modified  the  pressure  wave  generation  model  computer  urogram 
(Reference  13)  to  calculate  panel  deflections.  The  variable  image  method  (Reference  11)  for 
the  fluid/structural  reaction  was  employed,  and  Ankeney’s  structural  equations  (Reference 
12)  were  used  for  structural  response.  Computer  code  predictions  of  panel  deflections  were 
compared  with  a  single  exit  panel  shot  in  Reference  11.  Agreement  between  code  prediction 
and  experiment  was  quite  good.  The  computer  code  was  called  the  Unimodal  Hydraulic  Ram 
Structural  Response  (UHRSR)  program.  A  user’s  manual  wcs  published  by  Ankeney 
(Reference  14). 

The  UHRSR  code  was  modified  by  Wada,  Hsitz,  an<J  Jacobson  (Reference  15)  to 
include  a  capability  to  calculate  the  response  of  orthutropic  path's  to  hydraulic  ram  pressure 
loading.  Computer  calculations  of  wall  deflection  and  strain  were  compared  with  entry  panel 
experimental  data.  A  general  lack  of  agreement  was  obtained.  The  predictions  of  panel 
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strain  wars  particularly  bad.  The  most  probabls  causa  was  ths  fact  that  panel  deflections  did 
not  look  like  ths  assumed  cosine  curve. 

Herlin  and  Avery  (Reference  16)  merged  the  UHRSR  and  BR-1FC  codes  to  produce 
the  Hydraulic  Ram  Structural  Response  (HRSR)  code.  The  BR-1FC  code  was  an  extension  of 
the  BR-1  code  to  include  analysis  of  orthotropic  stnnUual  elements.  Comparisons  of  the 
HRSR  code  calculations  with  experiments  were  performed  by  Bradley  and  Avery  (Reference 
17)  and  Jacobeon,  Heitx,  and  Yamane  (Reference  18).  Reasonable  agreement  was  usually 
obtained.  Discrepancies  were  primarily  attributed  to  the  following: 

1.  The  finite  element  grid  is  usually  chosen  too  coarse  because  the  code  is  very 
expensive  to  run.  Large  pressure  gradients  near  projectile  exit  or  entrance  points  are  not 
adequately  represented. 

2.  Cavitation  at  the  fluid/structure  interface  is  not  accurately  modeled. 

Fuel  tank  failure  is  assumed  in  the  HRSR  code  whenever  the  material  strain  exceeds 
the  ultimate  strain.  Bradley  and  Avery  (Reference  17)  have  shown  that  damage  propagation 
from  ballistic  perforations  can  be  obtained  at  much  lower  strains  because  of  fracture 
mechanisms.  Ball  and  Fahrenkrog  (Referenco  19)  attempted  to  use  fracture  mechanics  to 
calculate  the  lengths  of  the  cracks  in  damaged  panels.  They  were  handicapped  since  their  use 
of  the  piston  theory  of  the  fuel/structural  interaction  gave  predicted  stresses  that  were  too 
low.  However,  they  were  encouraged  by  results  obtained  by  normalisation  of  the  stresses. 

Rosenberg  and  Bless  (Reference  20)  have  derived  an  empirical  scaling  model  of 
hydraulic  ram  failure  based  solely  on  fracture  mechanics  considerations. 

This  work  was  intended  to  extend  the  HRSR  code  to  include  multiple  supersonic 
prqjectiles.  The  HRSR  code,  however,  proved  to  be  very  difficult  to  run  and  inconvenient  to 
use.  The  UHRSR  code  was  available,  and  modifications  had  already  been  incorporated. 
Extensions  were  made  to  the  fluid/structural  interaction  model  to  include  cavitation,  and  the 
structural  model  was  generalised  to  widen  its  application.  The  resulting  computer  program  is 
called  the  ERAM  code.  This  code  will  be  documented  at  another  time. 

The  new  developments  in  the  hydraulic  ram  theory  are  presented  in  this  report. 
Experiments  were  performed  using  nontumbling  cylindrical  prqjectiles  to  give  repeatable 
date  for  exit  and  entrance  panel  shots  at  0-  and  45-degree  obliquity.  One  of  the  shots,  an 
entrance  shot  at  45  degrees,  is  analysed  extensively.  Predictions  of  fluid  pressure  and 
aluminum  panel  strains  are  compared  to  experimental  data.  Plots  of  the  hydraulic  ram 
variables  within  the  fluid  and  at  the  fluid/structure  interface  are  made  to  illustrate  the 
dynamic  interaction  of  the  fluid  mechanics  and  the  fuel  cell  structure. 

Fracture  mechanics  techniques  are  applied  to  predict  panel  catastrophic  failure 
thresholds  with  some  success. 
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PRESSURE  WAVE  GENERATION 

i 

Tht  theory  for  the  pressure  wow*  generation  for  dingle  tubeonic  penetratore  during 
the  drag  pheae  of  hydraulic  ram  wet  presented  earlier  in  References  3  and  4.  The  user’s 
manual  for  a  computer  code  implementing  the  theory  (Reference  13)  contains  the  derivation 
additional  equations  required  for  the  numerical  evaluation  of  hydraulic  ram  pressures.  The 
theory  has  been  extended  to  multiple  supersonic  projectiles.  The  intent  of  this  section  is  to 
document  the  extension  to  supersonic  projectiles  in  some  detail,  but  without  reproducing  all  of 
the  mathematics  for  subsonic  projectiles. 

The  hydraulic  ram  theory  assumes  that  the  flow  field  can  be  described  in  terms  of  a 
potential  function,  $,  which  satisfies  the  wave  equation 


csar 


(i) 


where  c  is  the  sound  speed  in  the  fluid.  In  terms  of  the  potential,  the  fluid  velocity,  n,  is  given 
by 


u  =  V<p 

and  the  pressure,  p,  is  obtained  from  Bernoulli’s  equation 


(2) 


1  J  (3 

P  =  P0-P—  ~  jpu 

where  Po  is  the  ambient  pressure  and  p  is  the  fluid  density. 

The  coordinate  system  for  the  analysis  is  shown  in  Figure  1.  The  coordinate  of  an 
arbitrary  point  is  (x,  w),  where  u  is  the  perpendicular  distance  to  the  trajectory  and  x  is  the 
distance  parallel  to  the  trajectory.  The  origin  of  the  coordinate  system  is  the  impact  point. 


FIGURE  1 .  Geometry  for  Pressure  Wave  Generation  Model. 
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The  simplifying  approach  of  tha  modal  is  to  try  to  approximata  tha  affact  of  the 
projectile  and  cavity  by  a  line  of  sources  distributed  along  the  traj  ectory  path  behind  the 
projectile.  Tha  potential  for  an  arbitrary  Une-source  distribution  is  given  by  the  integral 


f*4u 

♦(*.«,<)  *  - 


(4) 


where  t  is  the  distance  along  the  trajectory,  r  h  the  distance  from  0)  to  (x,  o>),  and  <  is  the 
source  strength  at  Since  4  satisfies  the  wave  equation,  the  source  strength  and  integration 
limits  must  be  evaluated  at  the  retarded  time  t  =  t- r/c. 


The  integral  in  Equation  4  is  a  line  integral  over  a  path  defined  by  i  =  t  -  r/c  with  the 
time,  t,  held  constant  Figure  2  shows  the  geometry  of  a  typical  integration  path  on  the  time- 
space  plane.  For  a  given  drag  force  history,  the  projectile  time  of  arrival  can  be  calculated  as  u 
function  of  projectile  position,  X^.  A  representative  trajectory  for  a  subsonic  penetrator  is 
included  in  Figure  2.  For  the  subsonic  case,  the  integration  path  intersects  the  trajectory  at  a 
maximum  of  one  point.  This  point  is  the  upper  integration  limit  X|»  which  is  the  projectile 
position  at  time  t  =  t^.  The  line-source  distribution  extends  to  the  impact  point,  so  the  lower 
integration  limit  is  X,  =  0. 

The  geometry  is  more  complicated  for  the  case  of  supersonic  penetrators.  The 
trajectory  and  integration  path  will  have  multiple  intersections  for  certain  choices  of  (x,  o>,  t). 


FIGURE  2.  Geometry  for  Typical  Integration  Path  on  Time-Space  Plane. 
Also  shown  is  the  trajectory  for  a  subsonic  penetrator. 
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This  situation  is  illustrated  in  Figure  3.  Three  integration  paths  are  shown  corresponding  to 
three  successive  times.  For  the  integration  path  labeled  t*  the  trajectory  is  tangent  to  the 
integration  path.  At  time  t*  pressure  disturbances  produced  by  the  projectile  first  reach  the 
point  (x,  <•>).  The  integration  limits  *  X*.  Because  of  the  tangency,  the  integral  will  be 
singular,  as  will  be  seen  below.  At  a  later  time,  ti  >  t,,  the  intersection  points  have 
separated.  One  integrates  between  limits  Xe(tfe)  >  X,(i*). 

For  subsonic  penetrators,  the  integral  (Equation  4)  for  the  fluid  potential  ia  taken 
between  limits  Xe(t)  and  0.  For  supersonic  penetrators,  the  integral  can  be  partitioned  into 
two  similar  integrals 


K  r 

Both  integrals  are  identical  to  the  subsonic  case,  so  that,  in  principle,  the  subsonic  results  are 
directly  applicable.  However,  required  derivatives  of  $  are  singular  at  the  sonic  point.  The 
singularities  arise  from  the  Doppler  shift  of  waves  produced  by  the  pcnetrator.  The 


CMSlANCe 


FIGURE  3.  Saquancs  of  Intagration  Curves  for  the  Potential  at  (x,  u)  for 
Times  t,  <  ti  <  tj.  The  trajectory  for  the  supersonic  panatrator  may 
intersect  the  intag raton  curve  twice,  as  shown  for  the  curve  ti. 
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singular! tiM  ara  rtaolvad  in  naiura  by  nonlinaar  terms  naglaetad  in  the  approximate  wave 
equation  for 

The  source  distribution  was  derived  in  References  S  and  13  as 


where  A(Q  is  the  maximum  cavity  radius  obtained  along  the  trajectory,  and  tb($)  is  the 
projectile  time  of  arrival  at  point  V  The  function  A($)  Is  given  by 

(-) 

AliJ“3wv 

where  the  derivative  (dfi/dXb)(  is  the  rate  at  which  the  projectile  kinetic  energy  is  deposited 
into  the  fluid  by  the  dreg  forces.  The  ambient  and  cavity  pressures  are  denoted  by  P0  and  Pc> 
respectively.  The  constant  B  has  dimensions  of  velocity  and  is  given  by 

»*  =  P°~Pc  (8) 
P  N 

where  N  is  a  geometric  dimersionless  parameter  discussed  in  Reference  3.  A  value  of  N  = 
1/0.434  is  used  in  the  computer  program. 

To  obtain  the  fluid  potential.  Equation  6  is  substituted  into  Equation  4  (or  Equation  5 
for  supersonic  penetrators),  which  could  then  be  integrated.  However,  for  the  desired  fluid 
pressure  and  velocity,  derivatives  of  the  potential  $  are  required.  In  References  3  and  13, 
expressions  for  these  derivatives  were  derived  as 


s*_  °A>  Li  *  +  \ 

at  ~21*.  V.x-JT  2  \x-Xl  +  RJ 


*1-^1- 
c  R 


a»  f  (  B 

~  V,  x-Xj  2  • 

c 


V 

a»  BAk  c  Rt  b 
ax  =  2 Rb  Vb  x-Xh  +  2 * 

1_  c 
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where 


r*»  « 

/„«* 

r* 

/,*  J^A-B<f-/fc)|i~cft 

For  supersonic  penetrators,  these  equation!  apply  to  derivative*  of  the  first  integral  in 
Bquation  5. 

The  integrate  I*»  and  Is  muat  be  evaluated  at  each  time  atep.  In  Reference  13,  fluid 
propertiee  were  calculated  at  a  sequence  of  time  etepe.  Hie  numerical  approach  wae  then  to 
calculate  time  derivativea  dlu/dt  and  4lx/dt,  which  were  then  numerically  integrated. 
Ai«alytical  expreeaiona  for  these  time  derivatives  were  derived  as 


Hus  approach  loses  accuracy  when  the  integration  time  step  is  large  or  the  penetrator 
velocity  is  dose  to  sonic.  A  more  satisfactory  approach  is  to  calculate  the  integrals  lw  and  I„  as 
functions  of  Xj,  instead  of  time  t  From  the  chain  rule  of  differentiation,  one  gets 


it  dXb  dt 

and  from  Reference  13  the  derivative 


(14) 


(15) 
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Solving  for  dl/dX^  and  substituting  Equations  12  or  13  and  Equation  IS,  one  gets 


Equations  16  and  17  are  numerically  integrated  to  obtain  Iu(Xt>)  and  Ix(X|,).  In  order  to  get 
Iu(t)  and  Ix(t)  required  for  Equations  10  and  1 1 ,  one  uses 


t  =  tb  +  -  \4x  -  Xb)2  +  to2 
In  Equations  9, 10,  and  1 1  a  singularity  occurs  when 


(18) 


V  -  6  -  c  <19 

*  R» 

As  shown  in  Figure  4,  the  left-hand  side  is  the  component  of  penetrator  velocity  along  a  line 
connecting  the  projectile  and  the  observer  at  (x,  w).  The  singularity  occurs  when  this  velocity 
is  sonic.  According  to  Equation  19  the  singularity  occurs  only  in  the  forward  hemisphere  from 
the  projectile.  In  the  computer  program  based  on  this  analysis,  the  singularity  was 
artificially  removed  by  imposing  the  limit 


1-t(tt)I  >0‘ 

When  this  analysis  was  initiated  for  supersonic  projectiles  it  was  thought  that  the  singularity 
arising  from  the  first  integral  in  Equation  5  might  be  canceled  by  the  second  integral  in 
Equation  5  because  of  their  opposite  signs.  This  turned  out  not  to  be  the  case. 

Pressure  waves  reflected  from  fuel  cell  walls  have  a  significant  effect  on  the  fluid 
pressure  distribution.  In  Reference  3  it  was  shown  that  the  typically  light  construction  of 
aircraft  fuel  cell  walls  produces  wave  reflections  as  if  the  walls  were  free  surfaces.  Since  the 
flow  is  described  by  a  potential,  the  wave  reflections  can  be  calculated  using  the  method  of 
images.  The  computer  program  described  in  Reference  13  uses  an  array  of  images  to  calculate 
the  wave  reflections  from  walls  of  rectangular  volumes.  The  surfaces  must  be  either  free 
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FIGURE  4.  Geometry  Showing 
Velocity  Components  for  Sonic 
Singulerity. 


surfaces  (using  negative  images)  or  rigid  surfaces  (using  positive  images).  The  surfaces  must 
all  be  of  one  type  or  the  other.  The  computer  program  based  on  this  analysis  has  a  slightly 
enhanced  capability.  Each  wall  may  be  individually  chosen  free  or  rigid. 


FLUID/STRUCTURE  INTERACTION 

The  analytical  prediction  of  hydraulic  ram  structural  effects  is  complicated  by  the  fact 
that  motion  of  the  structure  generates  a  pressure  within  the  fluid.  The  strength  of  the 
interaction  depends  on  the  fluid  density  and  is  very  important  in  hydraulic  ram,  whereas  it 
can  be  neglected  in  most  air  blast  problems  since  the  air  density  is  comparatively  low. 

The  fluid  pressure  acting  on  the  fluid/structure  interface  can  be  decomposed  into 
incident  and  reflected  pressure  waves.  The  incident  wave  arises  from  action  of  the  ballistic 
penetrator  and  may  include  wave  reflections  from  other  wall  surfaces  of  the  tank.  The 
incident  wave  is  independent  of  the  motion  of  the  fluid/structure  interface.  It  is  calculated 
using  equations  given  in  the  Pressure  Wave  Generation  section.  The  reflected  pressure  wave 
depends  both  on  the  incident  wave  and  the  motion  of  the  fluid/structure  interface.  To 
calculate  the  reflected  wave  involves  the  solution  of  a  diffraction  problem.  The  formal 
solution  is  known  and  requires  the  evaluation  of  KirkofPs  integral  (Reference  3).  Because  the 
approach  is  computationally  laborious,  several  approximate  solution  methods  have  been 
applied. 

The  variable  image  method  derived  in  Reference  11  is  used  here.  It  is  extended  to 
more  general  interface  boundary  conditions  so  that  more  complex  fuel  tank  wall  constructions 
can  be  examined. 
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Following  Reference  11,  the  pressure  Reid  in  the  fuel  cell  is  described  by  a  potential  <)> 
that  satisfies  the  wave  equation.  The  pressure  is  given  by  Bernoulli’s  equation 


P  = 


(21) 


and  the  fluid  particle  velocity  by 


u  ~  V<|> 


(22) 


The  fluid  density  is  denoted  by  p.  In  most  cases,  the  dynamic  pressure,  1/2  pu2,  in  Bernoulli’s 
equation  is  small  and  will  be  neglected  here. 


Let  the  potential  of  the  incident  wave  be  and  write  the  potential  of  the  reflected 

wave  as 


4>r  =  Q  ♦i 


(23) 


where  $i  is  the  mirror  image  of  $i  about  the  plane  of  the  wall.  In  general,  Q  varies  in  space 
and  time  such  that  $r  satisfies  the  wave  equation.  The  pressure  of  tne  reflected  wave  is  given 
by 


dQ 


P  =  -p —  =  —  pQ —  -  p«J>  — 

r  dt  at  KY'  at 

and  the  normal  component  of  the  fluid  velocity  is 


.  r  dQ 


u  =  —  =  Q  — -  +  <p  — 

rn  dn  an  ldn 


(24) 


(25) 


At  the  fluid/structure  interface  one  has 
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of  motion  of  the  wall.  Equations  31  and  32  are  used  for  the  fluid/structure  interaction  in  the 
HRSR  computer  program  (Reference  16). 

The  initial  value  of  Q  =  1  is  choren  in  order  for  Equation  31  to  give  the  correct 
pressure  doubling  effect  A  difficulty  with  Equation  32  arises  because  initially  $  =  o  and  uu 
=  o  and  is  therefore  indeterminant.  In  References  11  and  16,  this  indeterminacy  was 
removed  by  setting  the  initial  value  of  $  to  a  small  positive  value. 

A  second  difficulty  with  Equation  32  arises  when  the  wall  begins  to  move  before  the 
incident  wave  arrives.  Under  this  circumstance,  u«  o  and  <|>  =  o,  making  Equation  32 
singular.  In  Reference  16  the  approach  chosen  was  to  set  Q  =  1  before  the  incident  wave 
arrives,  then  use  Equation  32  thereafter. 

In  general,  the  fluid  does  not  remain  in  contact  with  the  fuel  cell  walls.  Large 
negative  pressures  are  often  developed  at  the  wall,  which  can  lead  to  fluid  cavitation. 
Equations  31  and  32  do  not  apply  in  this  situation.  The  approach  to  cavitation  taken  in 
Reference  16  was  simply  to  truncate  the  total  wall  pressure  at  zero. 

Cavitation  is  one  example  of  a  pressure-dependent  boundary  condition.  Another 
example  occurs  when  a  layer  of  crushable  foam  is  used  to  isolate  the  tank  structure  from 
hydraulic  ram  pressure  loading.  For  the  analysis  of  these  cases.  Equation  30  can  be  written 
as 


a  =  (1  -  Q)u  .  +  —1(1  +  Q)P.  -  PI 
n  '**  pc  1 

(33) 

dQ  1 

~ ■  =  — [(1  +<?)P  -P| 
dt  p$.  * 

(34) 

In  Equation  34,  P  may  be  a  constant,  as  in  the  case  of  cavitation,  or  may  be  a 
prescribed  function  of  foam  thickness.  Given  the  value  of  Q  by  integration  of  Equation  34,  the 
fluid  surface  velocity  is  then  calculated  from  Equation  33.  For  either  cavitation  or  foam 
layers,  the  distance  between  the  fluid  surface  and  the  wall  must  be  known.  The  fluid  surface 
position  is  calculated  by  integrating  Equation  33. 


STRUCTURAL  RESP0N3E  MODELS 


LAGRANGIAN  FORMULATION 

The  structural  response  models  are  derived  using  the  Lagrange  equations  of  motion. 
In  this  method,  the  panel  deflection  profile  and  strain  distribution  a-e  postulated  a  priori  as 
simple  functions  over  the  panel  surface.  The  functions  contain  one  or  more  parameters  that 
may  vary  with  time.  These  parameters  are  the  generalized  coordinates  of  the  Lagrangian 
theory  and  are  obtained  by  solving  the  Lagrangian  differential  equations  of  motion. 
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The  LiXgrangian,  L,  is  defined  in  terms  of  the  kinetic  energy,  T,  and  potential  energy, 
V,  of  the  panel  as 


L  =  T(q,q)-V(q)  (35 

where  q  is  a  vector  of  generalised  coordinates  having  components  q^.  The  vector,  q,  is  the 
time  derivative  of  q. 

The  Lagrangian  equations  of  motion  are  given  by 


d  f  dL  \  aL 
dt\dqk)~  3 qh  =Q* 


(36) 


or,  substituting  Equation  35  for  L, 


d/ar\_ar  av 

dt\dqh)~  dqk  dqt  + 

The  panel  motion  is  described  in  terms  of  the  deflection  v»,  normal  to  the  panel 


(37) 


w  »  u>(q,xh y)  v 

where  x  and  y  are  coordinates  spanning  the  panel,  in-plane  displacements  are  given  by  u  and 
v  in  x  and  y  directions, 


u  =  u(o,xj) 
v  =  u(q,xj) 

The  kinetic  ener?v  of  the  panel  is  given  by  the  integral 


(39) 

*0) 


T  = 


(41) 


where  y  is  the  mass  density  of  the  panel  material  and  h  is  the  panel  thickness.  The  integral 
extends  over  the  surface  of  the  panel.  The  kinetic  energy  T  the  in-plane  displacements  is 
neglected. 


The  generalized  force  acting  on  the  panel  is  given  by  the  integral 


Q*  =  I J  PdTdxdy  {42) 

where  P  is  the  hydrostatic  pressure  acting  normal  to  the  panel.  In-plane  components  of  the 
pressure  loading  acting  during  large  deflections  are  neglected. 
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where  the  bending  coefficient*  are  given  by 


E  'h* 

Sy'h* 

X 

12  ’ 

Dy=  nr 

G'h3 

(47) 

12  ’ 

D*y~  12 

Following  Reference  13,  the  bending  strain  energy  of  an  orthotropic  plate  element  is  given  by 


1  /  a *»»  £w  „  fw\1J 

Substituting  Equation  46  for  the  bending  moments  and  integrating  gives 

'  if  to  \2  i?w  3*u> 


fw  \ 


1  f  (9u)  \«  aw  Srw 

V  =  -  D  — -  )  +  2J&.— r  — 

fc  2  l  *\  a**  /  a*2  ay2 


d,{^J +^(£)V 

The  membrane  strain  energy  of  an  orthotropic  plate  element  is  given  by 


V  =  -h 
m  g 


eso,  +  V,  +  VJd*dy 


For  large  deflections,  the  strains  are  given 

du  1  (iw 
E*  ”  dr  +  2  V  a* 

d  v  1  /  dw  V 

ey“  ay  +  2Vc|y  / 


)’ 


du  do  dwdu; 

Y*y  “  ay  a*  at  ay 


Substituting  Equations  43  and  51  into  50  and  expanding  the  result  yields 


(48) 


(49) 


(50) 


(51) 
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V  -i 

H»  2 


+  B  *  a2  +  u  w*  +  7  u»4 
y  y  y  y  4  y 


+  £*  2u  u  +  u  w*  +  v  u»*  +  7U»au»8 
|  *  y  *  y  j*  2  *  *1 


+  G  I  u2  +  o2  +  2u  t*  +2u  UMti 
xr|  y  *  y  *  y  *  y 


+  2o 


w  ui  +  «*V  | 

1  »  j  *  yjj 


c txdy 


The  subscript  notation  is  used  for  the  partial  derivatives  of  u,  v,  and  w. 


(52) 


DISPLACEMENT  FUNCTIONS 

There  ere  several  possibilities  for  the  displacement  functions  (Equations  38,  39,  and 
40).  In  Reference  11,  the  equations  were  applied  to  an  isotropic  rectangular  panel  with 
dimensions  a  and  b.  The  displacement  functions  used  are 


w  =  w  cos  X  cos  Y 
o 

u  =  csin  2X  cos  Y 
o  =  ccosX  sin2Y 

where  the  arguments  X  and  Y  are  given  by 


(53) 


The  origin  of  x  and  y  are  in  the  center  of  the  panel.  The  amplitudes  w0  and  c  are  identified  as 
the  generalised  coordinates  q  in  Equations  38,  39,  and  40.  The  Lagrangian  Equation  37  is 
integrated  to  obtain  w0  and  c  as  a  function  of  time.  In  Reference  11,  the  structural  model  was 
restricted  to  isotropic  panels.  When  comparing  the  computed  w0(t)  to  the  results  of  a 
hydraulic  ram  exit  panel  test,  agreement  was  quite  reasonable. 


19 


NWCTP6770 


The  aquations  for  the  mors  gtnaral  cast  of  orthotropic  materials  were  developed  in 
Reference  IS.  The  displacement  functions  (Equation  S3)  were  used.  Predictions  of  panel 
stress  were  made  for  impact  walls  of  hydraulic  ram  tests.  Agreement  was  so  poor  that  the 
technique  was  abandoned  for  several  years. 

The  Lagrangian  equation*  described  here  give  a  satisfactory  approximation  to  the 
panel  motion  only  when  the  assumed  displacement  functions  resemble  the  actual  panel 
displacement  In  the  exit  wall  shot  of  Reference  11.  Equation  S3  was  adequate  to  show 
reasonable  agreement  with  wall  deflection.  However,  this  does  not  imply  that  the  calculation 
of  stress  will  give  similarly  good  agreement  The  shots  in  Reference  15  were  performed  on 
entrance  panels  that  are  expected  to  have  a  more  highly  concentrated  pressure  loading  than 
exit  wall  shots  and  therefore  a  different  wall  displacement  Also,  the  entrance  panels  in 
Reference  15  were  mounted  in  a  frame  that  allowed  in-plane  motion  at  the  panel  edges.  This 
boundary  condition  is  not  satisfied  by  Equation  S3,  which  assumes  clamped  edges. 

The  value  of  the  Lagrangian  technique  is  the  speedy  solution  of  the  structural 
response  equations  compared  to  a  finite  element  solution.  The  drawback  is  the  inaccuracy 
obtained  when  the  postulated  displacement  functions  do  not  adequately  approximate  the 
actual  panel  motion.  Three  different  sets  of  displacement  functions  are  included  here  so  thut 
one  may  choose  the  most  appropriate  for  a  given  structural  configuration  and  shotline 
geometry. 

The  displacement  functions  used  in  Reference  11  are  generalised  slightly  to 

iv  =  ui  cosXcoaY 

O 

u  —  u  sin  2X  cos  Y 

O 

o  =  v  cos X  sin  2Y 
o 

where,  again, 


The  amplitudes  Uo  and  v0  are  allowed  to  v\ry  independently  to  give,  hopefully,  a  more 
accurate  result.  Equation  54  should  he  ..sed  when  panels  have  fixed  boundaries  and  panel 
motion  resembling  the  first  mode  is  anticipated.  This  set  of  functions  will  he  referred  to  as  the 
fixed-boundary  displacements. 

The  test  .ixture  used  in  Reference  15  inspired  the  second  test  of  displacement 
functions 

io  =  u>  cos  X  cos  Y 
o 

u  -  u  sin  X  cos  Y  +  — ■w2sin2X  cosaV 
•  4a  0 


2U 
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u  -  v  cos  Y  tin X  +  -7  w*cos*X  sin  2Y 
•  46  • 

where,  again, 


r-3. 

o  6 

The  aecond  term  in  the  expressions  for  in-plane  displacements  (u  and  v)  is  included  to  make 
the  normal  streeaes  vanish  at  the  panel  edjgea.  Since  the  normal  stresses  vanish,  this  set  of 
displacement  functions  will  be  referred  to  as  the  (tee-boundary  displacements. 

Neither  the  fixed  nor  the  free  displacement  functions  are  appropriate  for  large  panels 
or  for  the  initial  stages  of  the  entrance  panel  problem.  Here,  a  Qaussian  set  of  displacement 
functions  is  defined  by 


where 


w  =  w 


♦  Y*) 


• = v + v+ ^ 


(56) 


X  ~ 


y-’o 

b 


The  constants  Su*.  Suy,  Sy*.  and  Sy?  determine  an  initial  uniform  stress  loading  applied  to  the 
panel.  The  Qaussian  amplitudes  (w0,  Uo,  vj,  as  well  as  the  Gaussian  widths  and  centers  (a,  b, 
Xo,  and  y0)  are  all  generalised  coordinates  satisfying  the  Lagrange  differential  equations. 


EVALUATION  OF  ENERGY  INTEGRALS 

The  integrals  for  kinetic  energy  and  strain  energy  need  to  be  evaluated.  This 
procedure  is  straightforward  but  very  tedious.  The  kinetic  energy  integral  (Equation  41) 
yields  the  same  result  for  both  cosine  law  displacement  functions  (Equations  54  and  55) 

T-  (57) 
8  0 

where  the  dot  denotes  differentiation  with  respect  to  time. 

The  Gaussian  displacement  functions  (Equation  56)  give  a  more  complicated  result 
since  the  generalized  coordinates  a,  b,  Xg,  and  y0  are  all  functions  of  time.  Substituting 
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Equation  56  into  the  kinetic  energy  integral  (Equation  41),  thon  integrating  ovar  an  infinite 
panal,  ona  gats 


nytafc  ( . «  3  3 

T  =  - I  ur  +  -  u>  —  +  -ur-r 

4  1  0  4  •»*  4 


if  *  •  l 

•  a  r*  .a  .6 

+  ie*—  +  ur-r  +  *•»**•  +  u> w - 
a*  6*  0  *  *6 

+  (58) 

S  *o6  J 

The  bending  strain  integral  (Equation  40)  also  gives  the  same  result  for  both  sets  of 
cosine  law  displacement  functions.  Substituting  either  Equation  54  or  Equation  55  into 
Equation  49  and  integrating  over  the  plate  area  gives  the  result 


_  jD  2D,  D  4D  \ 

*  16  A  ,*  oy  g*  0V  / 


Similarly,  substituting  the  Gaussian  displacement  functions  of  Equation  56  into  Equation  49 
and  integrating  over  an  infinite  plate  area  one  gets 


«_  /D  80.  D  40  \ 

v»=  +  Hv  +  i*  +  si5?) 


The  membrane  strain  energy  for  the  fixed  boundary,  cosine  displacement  function  set 
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The  Lagrange  equation  (Equation  37)  can  be  simplified  by  noting  that  T,  Vb,  and  Qk 
are  independent  of  the  in-plane  displacements  u0  and  v0.  For  the  case  of  Gaussian 
displacements,  they  are  also  independent  of  the  parameter  {).  Therefore,  for  these  variables, 
the  Lagrange  equations  reduce  to 


av  av  av 

m  hi  m 


du 


* 


=  0 


(64) 


The  membrane  energy  Vm  (Equationr  61,  62,  and  63)  for  the  three  displacement 
representations  can  be  written  in  the  quadratic  form 


V  =  C  +  C  u2  +  C  u  t>  +  C  v2 

m  o  uu  o  mi  o  u  uv  a 


+  C  u  u>2  +  C  v  w2 

UW  O  0  vw  o  o 


+  C  w2  +  C  wA 

W  O  UW  o 


(65) 


where  the  C’s  are  not  functions  of  the  Uq,  v0,  or  w0,  but  may  be  functions  of  p  for  the  Gaussian 
displacements. 

If  one  differentiates  Equation  65  according  to  Equation  64,  one  gets  the  system 

2C  u+Co+C  u/2  =  0 

UU  O  UV  O  UW  o 


The  solution  of  Equation  66  is 


C  u+2Cv+C  ur  =  0 

uv  o  wo  vw  o 


u  =  C  w 

o  u  o 


where 


v  =  C  w 

O  0  0 


C  =(C  C  -2 C  C  VA 

U  UU  VW  UV  UW 

C  =(C  C  -2C  C  yA 

v  UV  UW  uu  vw 


and  the  determinant,  A,  is  given  by 


(66) 


(68) 


A  =  4  C  C  -  C2 

uu  vv  uv 


(69) 


Thus,  Equation  67  gives  u0  and  v0  when  w,„  a,  and  b  are  known  (since  the  C's  are  functions  of 

a  and  b). 
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For  Gaussian  displacements,  the  Cu  and  Cv  are  functions  of  the  parameter  p.  At  this 
point  in  the  development,  the  last  of  Equation  64  is  still  not  satisfied.  The  variation  of 
membrane  energy  with  P  is  shown  in  Figure  5.  The  curves  in  Figure  5  result  from  Equation 
62  applied  to  an  isotropic  panel  with  a  =  b  and  Sux  =  Sy*  =  Suy  —  Syy  =  0.  Curves  are 
plotted  in  Figure  5  for  values  of  Poisson's  ratio  of  0.25  to  0.35,  which  covers  the  range  of 
interest  All  of  the  curves  show  a  broad  minimum  at  around  P  =  3.  This  value  is  fixed  for  all 
future  analysis.  The  other  minimum  in  Figure  5,  with  P  <  1,  corresponds  to  a  compressive 
membrane  stress  at  the  Gaussian  center  and  is  therefore  not  a  good  choice. 

The  distribution  of  membrane  strain  over  the  panel  is  given  in  Equation  51  as  a 
function  of  derivatives  of  displacements.  The  membrane  strains  are  calculated  and  plotted  in 
Figures  6  through  9  as  a  function  of  distance  along  the  x  axis  from  the  panel  center.  Along 
this  line  the  components  of  membrane  strain  eE  and  ey  in  Equation  51  are  labeled  radial  and 
hoop  strain  in  the  figures.  The  membrane  strains  are  normalized  for  an  isotropic  square 
panel  with  a  Poisson’s  ratio  of  0.3. 

Profiles  of  the  radial  and  hoop  components  of  the  normalized  membrane  strain  are 
plotted  in  Figures  6  and  7  for  the  Gaussian  displacements.  The  strains  are  plotted  for  several 
values  of  the  parameter  a  in  order  to  show  its  effect.  The  curves  corresponding  to  a  =  3,  which 
was  chosen  on  the  basis  of  minimum  energy,  appear  to  be  the  most  physically  realistic. 

Profiles  of  membrane  strain  for  cosine  deflections  are  plotted  in  Figures  8  and  9  for 
fixed  and  free  boundary  conditions,  respectively. 
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Exponent  for  inplane  Deflections 

FIGURE  5.  Normalized  Membrane  Energy  (Equation  63)  for  Gaussian  Displacements 
of  Isotropic  Square  Panels  as  «  Function  of  the  Parameter  p,  Which  is  Defined  by 
Equation  56. 
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Dimensionless  Radius,  R/a 

FIGURE  6.  Normalized  Membrane  Radial  Strain  Distribution  in  Isotropic  Square 
Panels  for  Gaussian  Deflections  Showing  the  Effect  of  the  Parameter  p. 


FIGURE  7.  Normalized  Membrane  Hoop  Strain  Distribution  in  Isotropic  Square 
Paneis  for  Gaussian  Deflections  Showing  the  Effect  of  the  Parameter  p. 
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Dimensionless  Distance,  2X/a 

FIGURE  8.  Normalized  Radial  and  Hoop  Membrane  Strain  for  Cosine  Deflections 
(Equation  54)  of  an  Isotropic  Square  Panel  With  Fixed  Boundary  Conditions. 
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EQUATIONS  OF  MOTION  FOR  COSINE  DEFLECTIONS 

For  the  cosine  deflection  functions  with  either  fixed  or  free  edges,  the  only  remaining 
generated  coordinate  is  w0.  To  evaluate  the  Lagrange  equation  (Equation  37),  first  rewrite 
the  kinetic  energy  integral  (Equation  57)  as 


where 


T=  ;Cu)S 
2  9  o 


r  _  V hab 
4 


(70) 


(71) 


The  various  derivatives  of  the  kinetic  energy  integral  required  for  Equation  37  are 


dw  p  0 

o 


±{Z.).CB 

dt\du>  /  p  • 


a t 

dW 


=  0 


(72) 


(73) 


Similarly,  the  bending  strain  energy  (Equation  59)  can  be  expressed  as 


V.  =  -  C.w 
b  2  6  0 


(74) 


where 


n  / D  2D,  D  4D  \ 


(75) 


whence,  the  required  derivative  is 


aV6 

=  C.w 
dw  b  0 

o 


(76) 
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The  derivative  of  the  membrane  strain  energy  is  obtained  from  Equation  65 


9w_ 


=  2C  u  w  +  2C  v  u?  +  4C  w* 

uw  o  o  umoo  mu  o 


(77) 


For  the  cosine  deflections,  the  Cw  term  vanishes,  so  it  was  not  included  in  Equation  77. 
The  Uf,  and  v0  in  Equation  77  can  be  eliminated  using  Equation  67,  whence 


-C  wa 

m  o 


(78) 


where 


C  =  2 C  C  +  2C  C  +4C 

m  uui  u  ihu  u  ww 


Substituting  the  derivatives  of  Equations  72, 73, 76.  and  78  into  Equation  78,  one  obtains 

C  w  +  C.w  +  C  =  Q 

p  o  oo  mo  w 


(79) 


(80) 


■> 

The  generalized  force  acting  on  the  panel  must  be  evaluated  numerically.  For  the  cosine 
deflections  it  is 

(81) 


Qu  =  j  J  PcoaXcoaY dxdy 

where  the  integral  is  taken  over  the  area  of  the  panel.  The  second-order  differential  equation 
for  w0  may  be  solved  by  a  variety  of  techniques  when  the  pressure  P  is  known. 


EQUATIONS  OF  MOTION  FOR  GAUSSIAN  DEFLECTIONS 

The  Lagrangian  equation  (Equation  37)  for  the  Gaussian  deflections  is  more  complex 
because  there  are  more  independent  generalized  coordinates.  The  in-plane  deflection 
equations  have  already  been  integrated,  resulting  in  Equation  67 .  There  remains  the  out-of¬ 
plane  deflection  (w0),  the  Gaussian  widths  (a  and  b),  and  the  Gaussian  centers  (x0  and  y0). 
Therefore,  five  Lagrangian  equations  are  required. 

The  total  time  derivative  in  Equation  37  is  expanded  in  the  usual  way  with  the  chain 
rule  for  differentiation 
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pt>  evity,  we  rewrite  Equation  84  in  the  form 


sK  ST  )  =  +  "tf5  +  "la*1  +  A1 

o 

dt  (  ad  )  “  Mnwo  +  wa®  +  +  A2 

df#r\ 

)SV#  +  V  +  V^3 

5(?)-JfA^4 

o 

dtrr\ 

•'o 


(84) 


(85) 
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Therefore,  the  Lagrengian  Equation  37  can  be  written  aa 


a T 


M'vt  +  Mja  +  hfjb  =  ^ 
31  o  32  33  36 


dV 

m 

dw 

a 

dV 

m 

da 

da 

** 

dV 

m 

a b 

~  db 

-  A,  +Q 
I  »* 
o 


*<4,5,  =  “A,  + 


Af«Br.ss-A3+^ 

'o 


(86) 


Vanishing  terms  were  not  included  in  Equation  86.  The  remaining  derivatives  are  easy  to 
derive,  but  are  lengthy  and  will  not  be  shown  here.  It  should  be  noted  that  the  derivatives  of 
the  membrane  potential  energy,  Vm,  are  evaluated  at  constant  Uo  and  v0.  These  variables  are 
then  eliminated  after  the  differentiation  by  using  Equation  67. 


The  equations  of  motion  (Equation  86)  have  the  form  of  a  mass  matrix  multiplying  an 
acceleration  on  the  left-hand  side  and  a  force  on  the  right-hand  side.  The  w0,  a,  and  b 
accelerations  are  coupled.  To  integrate  these  equations,  the  accelerations  must  be  solved  for 
by  inverting  the  mass  matrix.  This  is  done  numerically.  It  should  be  noted  that  the  mass 
matrix  becomes  singular  when  =  0,  the  initial  condition.  This  can  be  rationalized  by 
noting  that  the  shape  of  the  Gaussian  is  independent  of  the  widths  a  and  b  when  wQ  =  0;  it 
vanishes  everywhere. 


EXPERIMENT 

Hydraulic  ram  testing  was  performed  to  obtain  data  for  comparison  with  the  analysis. 
The  target  test  tank  was  a  2-foot  cube  filled  with  water.  Projectiles  wore  steel  right  cylinders; 
they  were  0.5  inch  diameter,  1.5  inch  long,  and  weighed  38  grams.  The  advantage  of  these 
projectiles  is  that  they  are  stable  during  fluid  penetration,  and  therefore  produce  a  very 
repeatable  hydraulic  ram  effect 

Entrance  and  exit  panels  were  constructed  of  0.063-inch-thick  7075-T6  aluminum. 
The  12-inch  square  panels  were  clamped  to  the  test  tank  using  a  frame  that  had  a  10-inch- 
aquare  inside  dimension.  The  effective  panel  size  is,  therefore,  a  10-inch  square.  Strain  gages 
were  applied  to  the  panels  at  points  located  1.5  and  3  inches  above  center.  The  projectile 
trajectory  aimpoints  were  at  the  panel  centers  for  each  shot.  The  gages  were  applied  to  both 
the  inside  and  outside  of  the  panels.  Many  gage  failures  were  experienced.  The  leads  to  the 
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gages  ware  not  (Irmly  attached  to  the  panels.  The  leads  pulled  loose  from  the  gages, 
particularly  during  shots  with  high  projectile  impact  velocities. 

Fluid  pressure  was  measured  by  up  to  six  piesoelectric  pressure  gages.  They  were 
typically  locaUri  1 1/2  inches  above  the  expected  shot  line  and  about  3  inches  apart.  The  peak 
pressures  obtained  ranged  flrom  30,000  pei  down  to  200  pei.  The  same  gages  were  used  for  all 
of  the  tests;  consequently,  the  low  pressure  results  are  somewhat  noisy. 

A  summary  of  the  shot  matrix  and  overall  test  results  is  shown  in  Table  1 .  Since  the 
target  panels  were  7075-T6  aluminum,  a  very  brittle  alloy,  only  two  target  damage  types 
were  obtained.  At  low  velocities,  a  single  perforation  was  obtained.  At  higher  velocities, 
cracks  originating  at  the  perforation  extended  to  the  panel  edges. 


TABLE  1.  Shot  Matrix  and  Results  Summary. 


Shot  no. 

Panel 

Angle,  deg 

Velocity,  ft/s 

Result 

234 

Entrance 

0 

1019 

—a 

23S 

3195 

236 

2050 

F 

237 

2550 

F 

238 

Entrance 

45 

2045 

F 

239 

9S0 

- 

240 

1563 

F 

241 

Exit 

45 

1205 

F 

242 

862 

- 

243 

Exit 

0 

1208 

F 

244 

893 

F 

245 

1 

766 

F 

•  Minimal  damage. 

4  Catastrophic  failure. 
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DETAILED  ANALYSIS  OF  SHOT  239 

The  hydraulic  ram  analysis  is  applied  to  shot  239,  a  low  velocity  oblique  entrance 
shot.  The  panel  suffered  minimal  damage  and  the  strain  gages  survived  the  test.  The 
geometry  of  the  test  tank,  entrance  panel,  and  the  cylinder  shot  line  are  shown  in  Figure  10. 
The  tank's  X,  Y,  Z  coordinate  system  is  also  shown  on  the  figure. 

Five  pressure  gages  were  used  in  the  shot.  They  were  placed  1  1/2  inches  above  the 
shot  line  and  spaced  3  inches  apart.  The  first  gage  was  3  inches  from  the  impact  point.  The 
five  pressure  gage  records  are  the  noisy  traces  in  Figures  11  through  15.  For  comparison, 
calculated  pressure  curves  are  plotted  on  the  same  figures.  The  analysis  assumes  a  drag 
coefficient  of  0.82.  Rigid  surface  reflections  are  assumed  for  the  massive  bottom  and  rear 
walls  of  the  test  tank,  and  a  free-surface  reflection  is  assumed  for  the  top.  Wave  reflections 
from  the  other  surfaces  are  ignored.  No  other  assumptions  are  required.  Comparison  of 
experiment  with  theory  is  quite  good. 
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FIGURE  11.  Prtssur*  Gagt  Rtcord  No.  1  for  Shot  239. 
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FIGURE  1 2.  Pressure  Gage  Record  No.  2  for  Shot  239. 


Time,  milliseconds 

FIGURE  13.  Pressure  Gege  Record  No.  3  for  Shot  239. 
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10  IS  20 
Y,  inches 

FIGURE  1C.  Pressure  Plot  at  0.1  Millisecond  After  Impact.  Shot  239. 


Fluid  pressure  end  velocity  are  calculated  over  the  horisontal  plane  Y  =  constant, 
which  contains  the  shot  line  (refer  to  Figure  10).  In  Figure  16,  the  pressure  is  plotted 
corresponding  to  a  time  of  0.1  millisecond  alter  impact.  The  bottom  curve  is  the  pressure  on 
the  entrance  panel  and  the  top  curve  is  the  pressure  at  the  rear  wall.  In  this  three- 
dimensional  plot,  deflection  of  a  curve  in  the  upward  direction  indicates  a  positive  pressure  at 
that  point  The  preasure  scale  on  the  left  gives  the  pressure  corresponding  to  a  given 
deflection  in  a  curve.  In  Figure  16,  the  highest  pressure  plotted  is  about  1.6  kpsi.* 

The  pressure  wave  discontinuity  arising  from  impact  can  be  seen  propagating  away 
from  the  point  of  impact.  The  fluid  cavity  behind  the  projectile  can  be  seen  plotted  as  a  region 
of  sero  absolute  pressure.  The  pressure  is  highest  in  the  vicinity  of  the  projectile.  Figures  17 
and  18  show  the  pressure  distribution  at  0.5  and  1.3  milliseconds,  respectively,  after  impact. 
The  region  of  high  pressure  moves  with  the  projectile,  and  the  cavity  increases  in  sise.  After 
the  projectile  exits  the  tank,  the  pressure  drops  to  a  comparatively  low  level.  In  Figures  16-18 
the  pressure  wave  reflections  were  not  calculated. 

Fluid  velocity  at  0.1  and  1.0  millisecond  after  impact  is  plotted  in  Figures  19  and  20. 
The  quantity  actually  plotted  is  peu  where  p  is  fluid  density,  c  is  fluid  soui\d  speed,  and  u  is 
fluid  velocity.  The  product  has  the  dimensions  of  pressure.  It  enters  into  the  variable  image 


*  kpsi  means  thousand  pounds  per  square  inch. 
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FIGURE  19.  Fluid  Velocity  at  0.1  Millisecond  After  Impact.  Shot  239. 
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(Equation  31)  in  the  same  role  as  pressure.  The  velocity  plots  have  two  differences  from  the 
pressure  plots:  the  peak  velocity  is  found  on  the  cavity  surface,  not  just  near  the  projectile, 
and  the  high  fluid  velocity  remains  even  after  the  projectile  exits  the  fluid.  The  magnitude  of 
the  peak  velocity,  pcu,  is  about  three  times  the  pressure  peak.  According  to  Equation  31 ,  the 
fluid  velocity  is  going  to  play  an  important  role  in  the  wall  pressure  loading. 

Structural  reaction  of  the  entrance  panel  was  calculated  using  the  assumptioi  of  a 
Gaussian  panel  deflection.  The  resulting  peak  deflection  is  plotted  in  Figure  21.  It  rises 
monotonically  to  0.7  inch  during  the  5.0-millisecond  period  analyzed.  The  motion  of  the 
center  of  the  Gaussian  deflection  is  plotted  in  Figure  22.  Both  the  x  and  y  coordinates  are 
initially  at  5  inches,  the  impact  point  at  the  center  of  the  10-inch  square  panel.  The  x 
coordinate  follows  the  motion  of  the  projectile  until,  after  about  2  milliseconds,  it  lies  outside 
of  the  entrance  panel.  Because  of  the  symmetry  of  the  panel  pressure  loading,  the  y 
coordinate  of  the  Gaussian  center  is  not  expected  to  move,  and  it  does  not. 

The  widths  of  the  Gaussian  profile  in  the  x  and  y  directions  are  plotted  in  Figure  23. 
They  both  grow  rapidly  until  at  1  millisecond  they  exceed  the  panel  boundaries.  It  is 
important  to  remember  that  the  Gaussian  deflection  structural  model  represents  an  infinite 
plate  that  is  loaded  only  in  the  area  of  the  actual  panel. 

In  the  Lagrangian  formulation  of  structural  response,  the  wall  loading  is  expressed  by 
generalized  forces  defined  by  Equation  42.  Therefore,  the  most  appropriate  measure  of  wall 
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Time,  milliseconds 

FIGURE  22.  Motion  of  Center  of  Gaussian  Panel  Deflection. 


FIGURE  23.  Widths  of  Gaussian  Profile  in  x  and  y  Directions. 
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loading  isa  weighted  average  over  the  surface  of  the  wall.  For  any  quantity,  f,  the  weighted 
average,  f,  is  defined  as 


-Ml 

where 

HI 

The  integrals  are  taken  over  the  wall  surface 
functions. 

Weighted  averages  of  the  incident  pressure  and  normal  fluid  velocity  are  plotted  in 
Figure  24  as  a  function  of  time.  For  velocity,  the  quantity  pcu,  which  has  the  dimensions  of 
pressure,  is  plotted.  In  Reference  8,  the  Mpistonn  theory  of  the  fluid/structure  interaction  was 
shown  to  give  poor  correlation  of  theory  and  experiment.  The  piston  theory  is  appropriate 
when  P  as  pcu  which,  according  to  Figure  24,  is  not  valid.  After  an  initial  transient,  it  is 
apparent  that  the  fluid  velocity  provides  the  dominant  wall  loading  mechanism. 

The  weighted  average  of  the  total  wall  pressure  loading  is  shown  in  Figure  25.  The 
magnitude  of  the  wall  pressure  is  low  compared  to  the  incident  pressure  and  fluid  velocity 
because  motion  of  the  wall  causes  pressure  relief  at  the  wall.  The  jagged  appearance  of  the 


FIGURE  24.  Weighted  Averages  cf  Incident  Pressure  end  Normal  Fluid  Velocity. 


fwdxdy 

wdxdy 

and  the  w  are  the  wall  normal  displacement 
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Time,  milliseconds 

FIGURE  25.  Weighted  Average  of  Total  Wall  Pressure  Loading. 

pressure  curve  is  caused  by  cavitation  collapse  at  the  wall  surface.  The  collapse  process 
actually  produces  a  continuous  pressure,  but  subdivision  of  the  wall  into  numerous 
computational  elements  separates  the  collapse  pressure  into  a  series  of  pressure  pulses. 

The  growth  and  collapse  of  cavitation  at  the  wall  is  shown  in  Figure  26.  The  curve 
labeled  pressure  is  the  fraction  of  the  wall  surface  area  that  is  covered  with  cavitation  and  the 
constant  pressure  boundary  condition  is  satisfied.  The  curved  labeled  velocity  is  the  fraction 
of  the  wall  surface  area  wetted  by  the  fluid  and  where  the  fluid  velocity  must  equal  the  wall 
velocity.  The  fluid  curve,  labeled  central  cavity,  is  the  area  fraction  of  the  cavity  following 
behind  the  projectile.  Initially,  the  entire  panel  is  wetted.  The  intense  pressure  loading  at 
the  penetration  site  pushes  the  wall  away  from  the  fluid  over  the  remainder  of  the  panel  until 
at  about  0.5  millisecond  the  wall  is  almost  entirely  cavitated.  Thereafter,  the  fluid  catches  up 
to  the  wall  and  the  region  of  cavitation  shrinks. 

When  using  Gaussian  deflection  profiles,  wall  cavitation  will  be  exaggerated 
whenever  the  width  of  the  Gaussian  is  commensurate  with  the  width  of  the  wall.  This  results 
from  ignoring  the  clamped  boundary  condition  at  the  panel  edges. 

The  area  fraction  of  the  panel  covered  by  the  central  cavity  is  significant.  At  5.0 
milliseconds,  Figure  26  shows  that  the  central  cavity  has  grown  to  cover  40%  of  the  entire 
entry  wall. 

The  average  value  of  the  variable  image  theory  Q  function  is  plotted  in  Figure  27. 
The  initial  value  is  Q  =  1.  The  wall  initially  reflects  pressure  waves  like  a  rigid  wall.  When 
the  wall  moves  outward  in  response  to  the  initial  pressure,  the  Q  drops  very  rapidly.  When 
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FIGURE  26.  Growth  and  Collapse  of  Cavitation  at  Wall. 


Time,  milliseconds 


FIGURE  27.  Average  Value  of  Variable  Image  Theory  Q  Function. 
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Q  =  -1,  the  wall  surface  reflects  pressure  waves  like  a  free  surface.  The  value  of  Q  seems  to 
stabilise  at  about  Q  =  -0.5,  which  may  be  an  artifact  of  the  averaging  for  this  particular 
problem. 

The  peak  wall  membrane  stresses  calculated  at  the  Gaussian  center  are  &hown  in 
Figure  28.  The  peak  stress  is  larger  in  the  X  direction,  which  results  from  a  smaller  Gaussian 
width  in  this  direction.  The  maximum  stress  la  about  45  kpsi,  which  is  significantly  less  than 
the  yield  stress  of  73  kpsi  for  the  material.  The  panel  did  not  fail,  but  it  had  a  barely 
noticeable  permanent  deformation  near  the  impact  perforation.  This  deformation  would  have 
resulted  from  additional  bending  stresses  not  shown  in  Figure  28.  Wall  strain  energy  is 
shown  in  Figure  29.  Most  of  the  energy  is  in  the  form  of  membrane  stresses,  particularly  at 
later  times.  The  wall  mr  "s  is  so  low  that  the  kinetic  energy  is  never  significant. 

Geometry  of  the  entrance  panel  for  shot  239  showing  strain  gage  locations  is  shown  in 
Figure  30.  Table  2  gives  the  data  for  each  gage. 

The  experimental  strain  gage  data  are  shown  compared  to  corresponding  calculated 
strains  in  Figures  31  through  36.  It  is  obvious  from  the  figures  that  the  simple  Gaussian 
structural  deflection  is  not  suitable  to  predict  the  details  of  the  strain  records.  In  fact,  some  of 
the  grosser  features,  such  as  the  large  compression  obtained  experimentally  for  the  Y  strains 
on  the  wet  side  gages,  were  not  predicted  by  theory.  The  Gaussian  deflection  profile 
predictions  are,  however,  much  better  than  similar  predictions  using  the  cosine  deflection 
profiles  for  impact  walls. 
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TABLE  2.  Strain  Gage  Location  and 
Orientation  for  Shot  239. 


Gag* 

no. 

Y  Offset, 
in. 

Strain 

direction 

Side 

1 

3 

Wet 

2 

3 

Wet 

3 

1.5 

Wet 

4 

1.5 

Wet 

S 

3 

VS:* 

Dry 

6 

3 

m 

Dry 

FIGURE  30.  Entrance  Panel  Geometry  Showing 
Location  of  Strain  Gages  for  Shot  239. 
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Time,  milliseconds 

FIGURE  31 .  Experimental  Strain  Gaga  Data  Compared  to  Corresponding 
Calculated  Strains.  No.  239*1. 


Time,  milliseconds 


FIGURE  32.  Experimental  Strain  Gage  Data  Compared  to  Corresponding 
Calculated  Strains.  No.  239*2. 
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FIGURE  33.  Experimental  Strain  Gage  Data  Compared  to  Corresponding 
Calculated  Strains.  No.  239-3. 


FIGURE  34.  Experimental  Strain  Gage  Data  Compared  to  Corresponding 
Calculated  Strains.  No.  239-4. 
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FIGURE  35.  Experimental  Strain  Gage  Data  Compared  to  Corresponding 
Calculated  Strains,  No.  239-5. 


Time,  milliseconds 


FIGURE  3€.  Experimental  Strain  Gage  Data  Compared  to  Corresponding 
Calculated  Strains.  No'.  239-6. 
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The  igrMnwnl  between  experiment  and  theory  it  the  beet  that  hae  been  obtained  to 
date  with  the  simple  structural  theories.  A  more  refined  calculation  using  a  finite  element 
model  would  explain  the  discrepancies  obtained  in  this  report 

In  order  to  understand  the  hydraulic  ram  phenomena,  the  special,  as  well  as 
temporal,  variation  of  the  fluid  and  wall  panel  properties  over  the  surface  of  the  wall  are  of 
interest  Plots  showing  special  distribution  are  shown  at  incident  pressure  (Figure  37), 
incident  normal  velocity  (Figure  38),  total  wall  pressure  (Figure  39),  panel  deflection 
(Figure  40),  wall  cavitation  (Figure  41),  and  the  variable  image  Q  function  (Figure  42). 
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(d)  Time  ■  0.858  millisecond 


FIGURE  37.  Spatial  Distribution  of  the  Incident  Pressure. 
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FIGURE  39.  Spatial  Distribution  of  the  Total  Wall  Pressure. 


53 


NWCTP6770 


12 

tit 


0  2  4  I  I  10  12 

Y,  indies 

(a)  Time  ■  0.0969  millisecond 


ti  1 1 «  i » ■  ■  1 1  1 1  i  ■  1 1 1  i  i.i  i  i  i.j  1 1  1 1 1 1  i  i 

0  2  4  6  0  10  12 


Y,  indies 

(c)  Time  ■  0.478  millisecond 


12 

in 


Y,  indies 


(b)  Time  ■  0.192  millisecond 


10 1  - - 

Ei  ■  *  ■  i  ■ » ■  ■  i » ■  ■  ■  i  >  t  ■ « i  ■  ■  ■  *~i  ■  » *  ■  i 

0  2  4  6  8  10  12 


Y,  inches 

(d)  Time  ■  0.858  millisecond 


FIGURE  40.  Spatial  Distribution  of  the  Panel  Deflection. 
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(c)  Time  -  0.478  millisecond  (d)  Time  *  0.858  millisecond 

FIGURE  41 .  Spatial  Distribution  of  the  Wall  Cavitation. 
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FIGURE  42.  Spatial  Distribution  of  the  Variable  Image  Q  Function. 


WALL  FAILURE  PREDICTION 

The  7075-T6  aluminum  alloy  used  in  the  experiments  is  a  very  brittle  material 
susceptible  to  crack  propagation.  In  the  hydraulic  ram  testing,  little  or  no  damage 
enhancement  is  obtained  at  low  projectile  impact  velocities.  As  the  velocity  is  increased,  a 
threshold  is  passed,  and  the  damage  increases  catastrophically  by  means  of  fractures  that 
propagate  from  the  perforation  to  the  edge  of  the  panel.  As  the  projectile  velocity  is 
increased  further,  the  number  of  propagating  cracks  increases. 


56 


■W*  ‘fQftm>nrBt  w*-- 


NWCTP6770 


A  slight  modification  of  classical  linear  elastic  fracture  mechanics  has  been 
successfully  used  to  predict  the  residual  structural  capability  of  ballistically  damaged  panel 
elements.  This  work  is  summarized  in  Reference  17.  Panels  are  ballistically  damaged  and 
then  stressed  to  failure  on  a  tensile  or  compression  testing  machine.  The  damage  is 
characterized  by  the  maximum  extent  of  the  damage  perpendicular  to  the  applied  stress.  This 
distance  is  termed  the  maximum  transverse  lateral  damage  (TLD).  The  critical  stress  for 
crack  propagation,  o,.,  is  given  by  the  expression 

is  (87) 

o  =  A  (nTLD/2r 172 

C  C 

The  parameter  Ac  is  an  effective  critical  stress  intensity  factor.  For  initial  damage  in  the 
form  of  a  sharp  crack,  A«  is  the  classical  stress  intensity  factor  Kc.  In  general,  Ac  is  slightly 
larger  than  Kc  and  depends  on  the  panel  thickness.  Reference  17  gives  details  and  data  for 
specific  material  and  penetrator  combinations. 

It  is  proposed  to  use  Equation  87  to  correlate  wall  panel  failures  for  hydraulic  ram 
where  Og  is  identified  with  the  peak  membrane  stress.  For  example,  for  shot  239  the  peak 
membrane  stress  is  about  45  kpsi  (see  Figure  28).  The  TLD  is  the  geometric  size  of  the  right 
cylindrical  projectile,  taking  into  account  the  effect  of  impact  obliquity  on  the  hole  size. 

To  evaluate  Equation  87  for  hydraulic  ram  damage  prediction,  the  peak  stress  was 
calculated  as  a  function  of  impact  velocity  for  the  four  engagement  geometries,  entrance  and 
exit  shots  at  0-  and  45-degree  obliquity  (see  Table  1).  The  calculation  was  performed  for  both 
Gaussian  and  fixed-boundary  cosine  panel  deflection  profiles.  The  results  are  shown  in 
Figures  43  through  46.  In  each  figure,  a  line  is  drawn  at  the  stress  level  predicted  by 
Equation  87  for  panel  failure.  On  each  curve  in  the  figures,  an  open  symbol  represents  a  shot 
at  that  velocity  where  the  test  panel  survived.  Closed  symbols  represent  catastrophic  damage 
propagation.  If  Equation  87  holds,  then  all  of  the  open  symbols  would  lie  below  the  line  of 
predicted  panel  failure  and  all  of  the  closed  symbols  would  lie  above  it. 

For  entrance  shots  (Figures  43  and  44),  it  is  not  clear  whether  the  cosine  or  Gaussian 
deflection  profiles  are  more  appropriate.  One  would  expect  the  Gaussian  curve  to  be  best  for 
small  panels.  From  the  analysis  of  shot  239,  however,  it  appears  that  the  10-inch  panel  is 
borderline  for  this  choice.  For  0-degree  obliquity,  the  cosine  predicted  stress  curve  seems  too 
high,  although  another  data  point  or  two  would  clarify  this.  The  Gaussian  predicted  stress  is 
too  low;  this  is  reasonable  since  the  clamped  plate  boundary  condition  would  raise  the  stress 
levels  in  the  panel.  At  45-degree  obliquity,  no  distinction  between  the  cosine  and  Gaussian 
curves  can  be  made  because  of  the  sparsity  of  data  near  the  panel  failure  point.  It  should  be 
noted  that  the  entrance  wound  for  the  980- (Vs  shot  had  a  crack  that  extended  0.6  inch  away 
from  the  edge  of  the  projectile  geometric  perforation.  This  would  tend  to  lower  the  predicted 
failure  stress.  In  general,  the  TLD  is  larger  than  the  geometric  perforation,  particularly  near 
the  ballistic  limit  of  the  wall.  See  Reference  17  for  details. 

Similar  curves  and  data  points  for  exit  panels  are  shown  in  Figures  45  and  46, 
respectively,  for  0-  and  45-degree  obliquity.  In  both  cases  the  Gaussian  curves  are  too  low, 
and  apparently  underestimate  the  stress  level  at  the  exit  hole.  The  cosine  curves  are  better, 
particularly  for  45-degree  obliquity. 
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FIGURE  43.  Results  of  Peek  Stress  Calculation  Performed  for  Both 
Gaussian  and  Fixed-Boundary  Cosine  Panel  Deflection  Profile 
Entrance  Shot.  0-Degree  Obliquity. 
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FIGURE  44.  Results  of  Peak  Stress  Calculation  Performed  for  Both 
Gaussian  and  Fixed-Boundary  Cosine  Panel  Deflection  Profile 
Entrance  Shot  45-Degree  Obliquity. 
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FIGURE  45.  Results  of  Peak  Stress  Calculation  Performed  for  Both 
Gaussian  and  Fixed-Boundary  Cosine  Panel  Defection  Profile  Exit 
Panel.  0- Degree  Obliquity. 
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FIGURE  46.  Results  of  Peak  Stress  Calculation  Performed  for  Both 
Gaussian  and  Fixad-Boundary  Cosina  Panel  Deflection  Profile  Exit 
Panel,  45-Degree  Obliquity. 
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At  O-degree  obliquity,  both  curvet  give  tt  much  higher  street  level  then  expected, 
based  on  fracture  mechanics.  As  a  possible  explanation  of  this,  the  two  exit  wounds  at  lower 
velocities  (766  and  893  ft/s)  have  a  different  appearance  than  for  the  remaining  shots  in  the 
test  series.  In  all  the  other  tests,  a  clean  puncture  was  obtained  in  the  shape  of  the  projectile 
with  some  possible  limited  spalling  or  cracking  around  the  edges-s  plugging  penetration. 
For  the  two  low-velocity  shots,  the  panel  was  only  sheared  around  half  of  the  circumference  of 
the  projectile,  and  then  the  hole  was  apparently  enlarged  by  tearing  one  or  more  large  flaps  of 
material.  The  exit  velocity  of  the  projectile  may  be  approaching  the  ballistic  limit  for  the 
panel  where  the  interaction  between  the  wall  and  penetrator  is  significant.  This  interaction 
is  not  accounted  for  in  this  analysis. 

At  46-degree  obliquity,  the  cosine  curve  gives  a  good  estimate  of  the  failure  threshold. 
The  low-velocity  shot  exit  wound  was  enlarged  by  limited  cracking  and  spallation.  This 
would  tend  to  lower  the  predicted  failure  threshold. 


SUMMARY  AND  CONCLUSIONS 

The  pressure  wave  generation  model  developed  in  Reference  13  has  been  extended  to 
multiple  supersonic  projectiles.  The  extension  to  multiple  projectiles  is  based  on  linear 
superposition  of  fluid  potentials.  The  predictions  for  multiple  projectiles  are  expected  to  be  as 
good  as  those  for  single  projectiles;  however,  no  other  tests  have  been  performed  to  verify  this. 

The  predictions  for  supersonic  projectiles  are  invalid  in  the  forward  hemisphere  of  the 
projectile  because  of  a  mathematical  singularity  obtained  there.  Predictions  in  the  /ear 
hemisphere  may  be  of  use  for  entrance  wall  pressure  loading.  Preliminary  results  of  analysis 
of  tests  with  multiple  supersonic  fragments  (not  reported  here)  support  this  idea. 

The  variable  image  model  of  the  fluid/structure  interaction  has  been  generalized  to 
the  case  of  pressure-dependent  boundary  conditions  at  the  fluid/structure  interface.  This 
boundary  condition  applies  in  several  important  situations,  including  cavitation  at  the  wall. 

Fuel  tank  hardening  concepts,  including  air  gaps  and  foam  isolation  layers,  can  be  analyzed 
with  pressure-dependent  boundary  conditions. 

The  structural  reaction  model  of  References  12  and  IS  has  been  generalized  using 
Lagrange  formalism  to  include  more  accurate  wall  deflection  shapes  and  boundary 
conditions.  A  variable- width  Gaussian  shape  was  found  to  give  acceptable  comparison  of 
theory  with  strain  gage  records  of  entrance  wal  1  testing.  * 

Nine  instrumented  hydraulic  ram  tests  were  performed  to  provide  pressure  and  strain 
data  for  comparison  with  analysis.  A  detailed  analysis  of  an  entrance  panel  shot  at  45-degree  > 

obliquity  was  presented  comparing  theory  and  experiment.  Plots  of  the  pressure-wave 
propagation  within  the  fluid  were  presented.  Additional  plots  of  the  fluid  and  structural 
properties  at  the  interface  showed  details  of  the  hydraulic  ram  phenomena.  The  importance 
of  the  fluid  cavitation  at  the  interface  is  emphasized. 

Predictions  are  made  of  the  threshold  for  catastrophic  failure  of  fuel  cell  walls  using 
the  extension  of  fracture  mechanics  of  Reference  17.  The  predictions  were  moderately 
successful  when  compared  with  experimental  data.  More  detailed  structural  analysis 
methods  should  improve  the  accuracy  of  the  predictions. 
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